Hybridization effect in the Falicov-Kimball model 
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The extended Falicov-Kimball model which contains hybridization between localized level and 
conduction band has been studied using the exact diagonalization method. We obtain the single- 
particle density of states and the valency from this method and analyze the behaviors of valence 
transition and activation gap. We find that the valence transition is continuous due to hybridization, 
while the activation gap closes discontinuously. Moreover, the onset of the valence transition and 
the closing point of the activation gap mismatch each other. Our result is well consistent with the 
recent experiment on SmBe. 



I. INTRODUCTION 

During last several decades, extensive studies have 
been done to understand two-band models such as An- 
derson model, Kondo model and Falicov-Kimball 
model(FKM) |^. These models are especially useful to 
explain many features of lanthanides and actinides since 
these materials have broad conduction haTid{s,d-band) 
and well localized valence band(/-&an(i) near Fermi level. 
These materials show some typical strong correlation ef- 
fects, for example large effective mass, non-integer va- 
lency, and small activation gap. 

We can categorize various two-band models into two 
groups according to their interband interactions. The 
first group includes hybridization between valence and 
conduction band. Anderson model, periodic Ander- 
son model, and Kondo model are well-known exam- 
ples emphasizing hybridization. Models taking inter- 
band Coulomb interaction into account fall into the sec- 
ond group. This interband Coulomb interaction is often 
called Falicov-Kimball term since they considered this 
interaction most crucial in understanding many features 
on lanthanides and actinides. FKM has been studied to 
explain valence transition and metal-nonmetal transition 
in non- integer valent rare-earth compounds |^-^ . 

SmBg and many other non-integer valent systems show 
phase transitions in valence and electric property when 
temperature or pressure varies. Increasing temperature 
makes d-electron more itinerant and external or inter- 
nal pressure reduces lattice constant and enhances crys- 
tal field splitting. These two effects broaden conduction 
band and lower the energy of d-electron. As a result, 
/-level crosses d-band and the system undergoes valence 
transition. Variation of these parameters also changes 
the configuration of localized /-electron and itinerant d- 
electron, which may causes metal-nonmetal transition. 

These transitions have been explained using the 
hybridization gap theory in the periodic Anderson 
model ll| . According to this theory, the valence tran- 
sition is due to the overlap between /- and d-band, and 
the appearance of semiconducting gap is the result of the 
hybridization between two bands. The width of the semi- 
conductor gap is proportional to the characteristic energy 



scale, the so-called Kondo temperature Tk, and the gap 
disappears continuously as temperature or pressure in- 
creases. In other words, the hybridization gap theory 
predicts the second order metal-nonmetal transition. 

A generalized two-band model, the so-called extended 
Falicov-Kimball model(EFKM) considering both inter- 
band Coulomb interaction and hybridization was pro- 
posed. Though much efforts have been devoted to 
this model for last two decades, no settled result on 
the nature of EFKM exists as far as we know. The 
most controversial point in EFKM is its valence tran- 
sition behavior. It is well known that FKM shows dis- 
continuous valence transition as pressure changes ^ . 
Some authors argued that this discontinuity still remains 
for small hybridization p^-jl7|. Others, however, have 
reported that there is no discontinuous transition in 
EFKM ||l2|,|l^. The only consensus is that the valence 
transition behavior in EFKM is so sensitive to the ap- 
proximation used that only approximation-free solution 
can clarify above controversy. 

In this paper, we study EFKM to show the valence 
transition behavior when hybridization exists and to elu- 
cidate metal-nonmetal transition which has been given 
by Cooley et al. |Q, ^ recently. They measured for that 
the activation gap of SmBg, which is a well-known non- 
integer valent material, closes discontinuously. This re- 
sult makes the hybridization gap theory suspicious 
since it predicts continuous gap closing. 

The transition behavior shown by Cooley et al. is strik- 
ingly similar to that of Mott-Hubbard transition. It is 
natural to suppose that the Coulomb interaction may 
play an important role in SmBg and other non-integer 
valent materials. The simplest model, which includes 
Coulomb interaction between localized level and conduc- 
tion band, may be the Fahcov-KimbaU model(FKM) [|. 

Recently Farkasovsky ||,^ has tried to explain the 
experimental results of Cooley et al. 0| using FKM. 
He showed using the exact diagonalization method that 
the Mott-Hubbard type of metal-nonmetal transition ap- 
pears in FKM when Coulomb interaction is small. This 
fact suggests that FKM may be considered as a possible 
microscopic model for non-integer valent materials such 
as SmBg. 
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Since the valencem, i.e. the expectation value of f- 
electron number operator in the ground state, must be 
an integer in FKM, non-integer valency is the result of 
statistical average of the mixtures of different integer va- 
lency. To have intermediate valency, the system must 
include the hybridization between the localized level and 
conduction band. 

We use the Lanczos exact diagonalization method 
to remove the artifacts of approximation in studying this 
model. Due to the limitation of the numerical method, 
we treat 1-dimensional chain of 8 lattice sites. For this 
system, we first get the DOS of FKM depending on en- 
ergy scale €/, and introduce the effective gap Ae// to 
explain valence transition in terms of the DOS. We find 
that only continuous valence transition is allowed when 
hybridization exists. From the DOS of EFKM, we find 
that the activation gap shows sudden disappearance at 
a particular e/ as in the case of FKM, while the valence 
changes continuously. Moreover, the onset of the valence 
transition and the closing of activation gap show clear 
mismatch. This is well consistent with the experiment 
on SmBe @, ||[. 

we introduce briefly the models and 
We present our results on D OS, valence transi- 

Section IV is 



In section 
method 

tion, and the activation gap in section |III 



devoted to the summary and conclusion. 



II. MODEL 

The FKM is a model explaining the properties of 2- 
level system competing among binding of /-electron (e/), 
itinerancy of (i-electron(f), and on-site Coulomb interac- 
tion between /- and d-e\ection{U) . Ignoring spin degree 
of freedom, the Hamiltonian of FKM is written as 



Hn = -t 



(1) 



where d\{di) creates(annihilates) a spinless d-electron at 
i-site. n{ and nf denotes the i-site occupation number 
operator of /- and d-electron, respectively. 

In Eq. (|l|) , we must note that — ^ ^ n{ and n'^ = 
nf commute with the Hamiltonian. This commu- 
tativity implies that their expectation values, (n^) and 
(n''), are good quantum numbers of this model. There- 
fore {n^) and {n'^) remain some natural numbers from 
the beginning to the end, and all (nf ) and {nf) are or 
1. 

Due to this simplicity, there are several rigorous theo- 
rems and many numerical results about FKM (see ref. ||] 
and references therein). These are worthy of note since 
FKM can be applied to various systems having /- and d- 
electrons. For example, FKM has been regarded as a sim- 
plified version of the single band Hubbard model [^-p2[ , 
the crystalline formation model |^,|2^ and the binary al- 
loying model [p3[. 



Ignoring the hybridization, however, could lead to an 
oversimplification in describing non-integer valent sys- 
tems. In particular, it seems unreasonable to disregard 
hybridization in non-integral valent systems where the 
energy difference between /- and d- level is very small. 
If we remind that hybridization plays a key role in the 
periodic Anderson model, it is natural to add the hy- 
bridization term to Eq. (Q) . This leads us to an extended 
model 



vY,{fld, + d\h) 



(2) 



where V is the transition amplitude between /- and d- 
level. This is the Hamiltonian of EFKM. 

Using Hartree-Fork approximation, Avignon and 
Ghatak |l4j insisted that EFKM shows discontinuous 
valence transition as well as continuous one by decreas- 
ing the energy gap between two levels for a finite f-d 
hybridization. Da Silva and Falicov [|l5| also reported 
that they found first and second order valence transi- 
tions in EFKM using Green's function method with mean 
field approximation and zero conduction band width. 
According to their result, discontinuous valence tran- 
sition appears when U /V > 2.5. Kanda et al. |^] 
showed the same result with Da Silva and Falicov |1^] 
within Hartree-Fork approximation but non-zero conduc- 
tion band width. Recently, Giesekus and Falicov jl^ also 
reported using Hartree-Fork approximation that EFKM 
has first order valence transition for particular parameter 
range. 

However, there are some contradictory results to those 
mentioned above. Using Hartree-Fork approximation 
and taking the periodicity of the system into account, 
Leder p2| argued that there is no discontinuity in va- 
lence transition in EFKM. Baeck and CzychoU also 
reported using different decoupling method that first or- 
der valence transition does not exist in EFKM. 

Depending on the approximation used, some of results 
show discontinuity in valence transition but others do 
show only continuous transition. This is why we use the 
Lanczos exact diagonalization method to get rid of the 
effect of approximation. At first, we use the modified 
Lanczos method to get the ground state |2^]. Then us- 
ing the Lanczos method, we investigate the ground state 
valency {^pg\nf\^pg), and the DOS 

Z?oM = E \{^:!-'\%J^^)\'6{u;-E^+E^-') 

k,n 

= E^o(^,^) (3) 



and 



= J2S,ik,iu), (4) 
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where subscript 'o' and 'e' mean 'occupied' and 'empty', 
respectively. S{k,uj), the spectral function, means the 
partial DOS of momentum k. S{k,u!) can be compared 
with the angle resolved photoemission spectroscopy data. 



III. RESULTS 

The /- and c?-electron DOS of FKM depending on Cf 
are shown in Fig.|l|. We put in the figure the activation 
gap Aact corresponding to the energy gap between filled 
/-band and empty d-band and the effective gap A^ff 
corresponding to the energy gap between empty /-band 
and empty d-band. 
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FIG. 1. The DOS for FKM with V=0 and U=1.5. Fermi 
surfaces are located in between two large /-peaks in all cases. 

The ground state electron configuration in real space 
for 6/ < 0.47 is the alternating type such as -f-d-f-d-. 
There is no sign of change in electron configuration from 
ej=0(a), to e/=0. 47(c). The only change in DOS is shift- 
ing of /-DOS, to high-frequency side. On the contrary, 
Fig.|l| (d) is quite different from Fig.|] (a)-(c), which im- 
plies that sudden change takes place in between e/=:0.47 
and 6/ =0.48. This is clear from Fig.|| which shows /- 
electron valency (n/) depending on e/. In Figj^ (b) for 
[7=1.5, valence transition occurs at 0.47 < e/ < 0.48, 
where Ae//=0. 

Therefore it is clear that Ag/ / is the effective gap gov- 
erning the valence transition in FKM. The difference be- 
tween Aact and Ae// is resulted from onsite electron- 
hole binding and proportional to the onsite interband 
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FIG. 2. (n/) depending on for [7=0. 6(a), 1.5(b), and 
3.0(c). Hybridization effect is shown in each case with V=Q, 
0.02, and 0.06. Hybridization effect is more clear within large 
U or large e/ region. 



Coulomb interaction strength U . This interpretation is 
well consistent with the original idea of Falicov and Kim- 
ball 1^] about the origin of discontinuous valence transi- 
tion in two band system. Thus Ae// corresponds to the 
shrinking gap 



Aeff = Aact - 2(7(1 - (n/)) 



(5) 



suggested by Falicov and Kimball 

Fig.|| shows the type of valence transition. For V—Q 
the /-electron valency shows steps as £/ increases, which 
is the well-known discontinuous valence transition in 
FKM. The hybridization effect in the valence transition is 
also clearly shown in Fig.||. Within all parameter range, 
all edges of V=Q disappear for any small hybridization. 
This disappearance clarifies continuous valence transition 
in all U for V ^ {). It is interesting to see that the hy- 
bridization effect is more prominent in large U and large 
e / . This is consistent with the reduction of Ag/ / , which 
is corresponds to the length of plateau in Fig.0, as U or 
£/ increases. Eq.(H) clearly explains these behaviors of 

Ae//. 

The continuous transition can be confirmed via ana- 
lyzing DOS. We present the change of DOS with e/ for 
a finite hybridization V^=0.02 in Fig.|[ Comparing Fig.^ 
with Figj^, we find that small /-peak indicated by arrow 
in (c) grows as we increase e/ from to 0.47. This new 
creation of /-peak make Ae//=0. Therefore Ae// is not 
defined in the hybridization-induced valence transition 
or Ae// is always zero for F 7^ 0, and this means the 
transition is continuous. As this /-peak grows, valence 
electrons can transfer more easily to conduction band and 
we observe continuous change of valency as a result. 

We now study the transport property of the system in 
terms of the activation gap A^ct which has been mea- 
sured by Cooley et al. 0]. We calculate the activation 
gap Aact from DOS. We find that Aact decreases linearly 
as e / increases in the constant valence regions which cor- 
respond to plateaus in Fig.y. This is quite natural since 
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FIG. 4. Results on (a) activation gap Aact and (b) con- 
duction electron density (jii) depending on e/. ^=0.02 and 
?7=0.6. Line of (a) is drawn for eye guide. Our results show 
the sudden closing of Aact and the mismatching between the 
(jidj increasing and A^ct closing. These are well consistent 
with recent experimental results. 



FIG. 3. DOS for [7=1.5 and finite hybridization V=0.02. 
Fermi surfaces are located between two large /-peaks in all 
cases. Arrow in (c) indicates prominent difference between 
1^=0.02 and V=^ case which is shown in Fig.nl. 



the effect of e/ on DOS is shifting of the /-DOS with 
d-DOS as rigid background. Calculation has been done 
for 17=0.6 and y=0.02 and the result is shown in Fig.^. 
?7=0.6 is the value for which Farkasovsky has shown 
metal-nonmetal transition for V—^S |^ , 

Our result demonstrates that the activation gap closes 
not at the onset of valence transition but at the end of 
valence transition. This is the well-known mismatching 
showing in the experiment of Cooley et al. 



IV. CONCLUSION 



competition between Ag// and V. For a fixed V, the 
hybridization effect becomes more prominent as U or e/ 
increases. Figj^ clearly shows the increasing of J7 or e/ 
causes smoother transition. The DOS of Fig.|| shows 
that hybridization makes the valence transition contin- 
uous by showing the growing of small /-peak as e/ in- 
creases, which makes the effective gap vanishes. Con- 
tinuous growth of this /-peak reflects continuous valence 
transition when V ^ 0. 

Finally, we show that the exotic metal-nonmetal tran- 
sition behavior of SmBg [Q can be explained by EFKM. 
Due to finite hybridization the valence transition takes 
place continuously, while Aact decreases linearly as e/ in- 
creases until the system meets new valence phase. Aact 
closes suddenly at the point of valence change. These fea- 
tures of EFKM are well consistent with Cooley et al.'s 
experimental result on SmBg. 



In this paper, we consider a generalized two band 
model, the so-called extended Falicov-Kimball model 
which includes both the Coulomb repulsion and hy- 
bridization between valence band and conduction band. 
The behaviors of valence transition and metal-nonmetal 
transition in EFKM have been studied using the Lanczos 
method. 

First, we give a remark on the energy gaps. The activa- 
tion gap Aact governs transport properties of the system, 
while the effective gap A^ff governs the valence transi- 
tion. The latter is the effective gap suggested by Falicov 
and Kimball 1^ and determines the condition of valence 
transition in FKM. Intrasite electron-hole pairing causes 
decreasing of A^ct — Ag//. This is known as gap shrinking 
by Falicov and Kimball 

Secondly, we find that the valence transition becomes 
continuous even for small hybridization. The effect of 
hybridization in valence transition is determined by the 
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